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ABSTRACT 

Shape memory alloys (SMAs) received significant attention by biomedical, aerospace and automotive industries 
due to their captivating properti.es called shape memory effect (SME) and Pseudo elasticity. However, the applications of 
these alloys are mostly found to be biocompatible, and are very sensitive to aerospace and automotive applications. 
This realization of developing the materials amenable to both aerospace and automotive apptications along with 
biomedical apptications needs an assignment of developing suitable constitutive model, to understand the behavior of 
these materials in different environments. As a result, the present paper emphasizes a short review on various constitutive 
models, cited in the global scientific community allowing to understand the significance of shape memory alloys, and also 
suggests essential facts required for developing new models. 
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INTRODUCTION 

Recently, SMAs received huge attention by aerospace and automotive industries because of its intriguing 
properties, say Pseudo elasticity and SME. But, most of these SMAs are constrained to biomedical applications 
because of their compatibility; hence they cannot be directly integrated with aerospace and automotive industries. 
These industries require properties like high strength to weight ratio, high thermal resistance and corrosion 
resistance which requires different processing methodologies unlike those for biomedical industries. 
Present challenges in processing these alloys received a substantial attention because of unavailability of proper 
constitutive mechanics to identify the material behavior for different alloying composition. Nonetheless, significant 
work has been done for various processing methodologies influencing the properties of Nitinol SMA and 
developed constitutive relations in these aspects. The present section will emphasize few research works performed 
over the years in deduce the mechanism of these alloys, and also highlights the research gaps identified from the 
literature review. 

REVIEW OF CONSTITUTIVE MODELS 

During 1980’s, which is the early phases of articulating and modeling of the SME, engineers and 
scientists found a close analogy with ferroelectric bodies due to the similar hysteretic behavior. As a result, 
Miiller and Wilmanski (1) proposed a macroscopic model to understand the behavior of the hysteresis with 
ferroelectric bodies. The primary objective of Miiller and Wilmanski model is to simulate a body under uniaxial 
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tensile or compressive loading, by statistical mechanical approach to calculate the free energy and entropy of the system 
given by equation(l-l. Alongside, they also explained load-deformation and deformation-temperature curves, but the paper 
couldn’t explain the elastic modulus (E) at the origin against the temperature, as shown in Figure 1. However, the paper 
credits the modelling of lattice structure using statistical ensemble average model, which intended the basis for signifying 
the phase transitions in Pseudo elastic bodies. 

<P P ottD,T) = -NkT[ln2z) +flJ^-+ln 1 1 S 1 Q (l - coshpAdA ], 

<Pvot(P,T) = (E 0vot + Z^ = ^0(A)1V 0a .) - TH 0 (1-1) 





Figure 1: Load Deformation Curve for Intermediate Temperatures 
(C, D) and High Temperatures (E, F) 


Nevertheless, to ascertain the concept of free energy, an enhanced micro model is earmarked already on the basis 
of stress-strain curves and latent heat of phase transition with shear strain E by (2), which is well known for Falk’s Model. 
The only drawback Falk’s model account for is, constraining the martensite variants during phase transition by Schmid’s 
law (3). The Helmholtz free energy F per volume as a function of shear strain E and temperature T in Gibbs free energy is 
specified in 


G = F — X E ; rescaled to g = f — ae 

Where F{E, T) = ccE 6 - fE 4 + (ST - y)E 2 + F 0 (T ) 



Where oc, /? and y need to be determined experimentally, and at the same time, these values are completely 
dependent on the composition strictly. 


Jurgen Sprekels (4) proposed a similar model (2)for a one dimensional constitutive model using non-convex 
potentials of (5), which particularly connects with B2-M phase transitions. The free energy F is defined in 
Ginzburg-Landau form as given in equation (1-3. 
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F(e, e x , 9) = -C v 9 log (£) + C v 6 + C + k 1 (0- 6 ± )e 2 - k 2 e 4 + k 3 e 6 + r -e 2 (1-3) 

Where fcj ,k 2 ,k 3 , Y are constants to be derived, experimentally. From the drawbacks of Falk’s model, a rheological 
model of stress-temperature is constructed by (6) as shown in Figure 2, which accounts the equivalent direction and 
magnitude of critical lower shear stress from the energy equations. 


t^rf f—'gr—l 
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Figure 2: Rheological Model 

Where, material is entirely idealized by elastic springs (Hookean elasticity) and dry friction element (plastic/St.- 
Venant) defining temperature dependence stiffncss C(T), K(T ) and friction coefficients IJ-(T) given by relations in 
equation (1-4 


a Fl = K(T) + fi(T),a F2 = K(T) - p(T) 
Energy Stored in elastic spring -Csf t and in 


Pre-strained springs Ks pl (1-4) 

The only drawback that paper accounts for is the absence of strain hardening resulting in the following cases 


For low temperature K = 0 T in the absence of 

For high temperature p = 0 f strain hardening 


It is intriguing to understand intuitive nature of SMA to detect their yield point by simple stretch test. 
Theories like Gibbs free energy, Falk’s model and other models couldn’t explain the phenomena of yielding. 
General assumption of yield point from Gibbs energy concept suggests that discontinuities of deformation or stretch are 
observed, but the case of SMA doesn’t obey this concept because of its Pseudo elasticity. In this context, (7) proposed a 
model to predict the yield point by using Maxwell line with stress as a function of temperature given by equation 

d-5 

t* = r*(6 e ) (1-5) 

Later the work from (7), Falk (8) realized that domain wall motion is responsible for SME followed by the phase 
transformation, and hence investigated the domain wall motion under the influence of external shear forces. In his 
subsequent paper (9), he realized the importance of domain walls (where strain jumps) by assuming them as shock waves 
because of their small thickness. By introducing the shock waves into constitutive equation(l-2, shear stress is determined 
and is given by the equation (1-6 

a(E, T) = = 6 E 5 - 4 E 3 + 2y(T)E, a(E, S) = = 6 E 5 - 4 E 3 + 2y(S)E (1-6) 

oE oE 

It is important to know the magnitude of the phase transition, through which continuous-cooling-transformation 
(CCT) and time-temperature-transformation (TTT) diagram can be derived, a work in this pretext is already available by 
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constructing a constitutive model using continuum mechanics which can be reviewed from (10). This is the first 
macroscopic Phenomenological model developed to understand solid-solid phase transition and named as Tanaka model. 
The general thermo-constitutive equation for any material in a stress-strain state is given by 

_ = dc, jjff ^ = g f ^ 

d£ dT Sf ' ’ 


n d 0 

D =Po^ 

Whe re^0 = Po£ 


d 2 ^ 

il — Pn - 

F0 dEdl, 


(1-7) 


Substituting the internal energy in first law of thermodynamics, will result in equation 
(1-77, which yield the expression for rate of entropy as 



30\ 

s) e 


-(■ 




30 

i? 


^-Vv ur “> 0 

p 0 l 


-Po dj_ 

p dX 


( 1 - 8 ) 


The generalized form of equation 
known as Hoffman’s-Sprekel’s model (11). 

0 (£, O ) = 0 O ( 0 ) + 0 X ( 0)£ 2 + 0 2 ( e ) 


(1-8 is given below and is 

(1-9) 


Linear 1-D models have been developed so far to understand the SME, but during unloading at B2 state, doesn’t 
give the same path in hysteresis, but still behave like a complimentary shear. It indicates that there is energy dissipation 
which doesn’t allow the system amenable to apply basic second law of thermodynamics to understand SME. Therefore, a 
new term is added as dissipation state function to this law, and developed a non-equilibrium thermostatic model for 
stress-strain B2M path as given by 

Fbam ~F = (F BAMx - F x )exp (M0) 

I f BAMx~ f x I 

It is an imperative model developed by (12) which decomposes the strain tensor into elastic strain efj (obey Hooks 
law for B2 &M) and inelastic strain e“(govcrned by normality theory of plastic theory) given as 

Elastic e u = ef,- + e u , Inelastic de u = dA— (1-11) 

l J lJ IJ’ 1 J dSi i ' 



Extending the 

work of 

Tanaka 

model 

(10), Liang 

(13) expressed 

equation 

(1-7 

in integral 

form by 

assuming 

D,0,n 

as constants 

resulting in 

equation 


( 1 - 12 . 

= D(e- £q) + 0(T - T 0 ) + 22 (f - f 0 ) d-12) 


Unlike Tanka model, Liang considered Martensitic Volume fraction as an expression of cosine and exponential 
form, which enable the model to understand SME in SMA more clear form. The model is named after his name as Liang’s 
model. The above research is produced over the decade 1980-1990 is basically on the modelling of martensitic 
transformation temperature, shear stress, shear strain relations followed by slight emphasis on hysteretic behavior. But few 
models lack traces of experimental evidences, which were produced over the next decade i. e., 1990-2000. To discuss 
works on different parameters influencing the constitutive model, (14) established a critical experimental work to 
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understand factors effecting SME such as ageing of martensite on shape-strain degradation by subjecting to number of 
thermal cyclic loadings, but the work didn’t emphasize the mathematical modelling, and hence didn’t receive much 
attention. 

Otsuka (15) envisaged about the limits of stress and strain where complete strain/shape recovery can be attained. 
He explained this in terms of stress induced transformations (crystallography of martensitic transition) and its 
thermodynamics. The shape strain than can be written as P ± = RP 2 B , when resolved into parallel and normal to invariant 
plane such that, 

+ m^p^ (1-13) 

Where,m^ d p and ?n"pj represent shear and dilatational components respectively and the meaning of the equation 
(1-13. Since it is a plane strain problem, by using Mohr’s circle, the external work done can be expressed in terms of plane 
stress and habit plane orientation as 

A G s = -\o a \[m^sinZcp ± nt"(l + cos2<p)] (1-14) 

The thermodynamic stress-induced transformation than can be given as equation 
(1-15 called as Clausius-Clapeyron relation. 

3cr A S P ~ M _ A H* P ~ M _ AQ((7) 

dT ~ _ s p ~ M ~ ~ T 0 {o)£ p ~ m ~ ~ T 0 (p)e p ~ m 

Transformation strain is calculated based on the orientation dependence given as 

1 

£ = [( m^sinx ) 2 + 2mi!^sinxcosA + l]z — 1 + rn^sin 2 / (1-16) 

Waymen (12) developed multi-dimensional constitutive model using plastic flow theory whereas Liang (13) 
developed one dimensional constitutive model; using these two models a multi-dimensional constitutive model based on 
martensite fraction is developed and tested for nitinol SMA by (16). Using equation 
(1-7, the constitutive equation can be described by and validated with a SMA specimen subjected to torsion. A more 
generalized one dimensional constitutive model is developed by Spies (17) with inclusion of thermodynamic potential into 
Falk’s model (18) via variables, such as time-dependent distributed and boundary inputs, viscosity effects, fading thermal 
memory, local curvature effects and internal variables. This setup allowed the present model to solve all the parameters of 
Falk’s model. The well-posedness of the equation (1-3 is given in the dissertation of (17). The physical essence to the well 
posedness solution is gives by establishing a continuum model developed by Tanaka (10) and Liang (13) to Spies model 


(17). Inclusion of martensite internal variable 

(1-7, Brinson (19) proposed a new constitutive relation as follows 

function 

in 

equation 

{ds} = {ds^} + {ds*-} + {di p } + {ds T } 



(1-17) 

o — o 0 = D{e — £ 0 ) + o(r — r 0 ) + 12 s (f s — f s0 ) + n , T ( f T — <f 7-0 ) 



(1-18) 


A model developed by (20)which is inevitably more complicated form, which introduces phase change hardening 
that can examine partial transformation and transformation hysteresis. The journal follows the constitutive rate equation for 
elastic response in elastic state and phase transformation followed by the transformation laws. The martensite volume 
fractions expressed in Liang model (13) defined by equation 
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(1-7 is derived here as 

—T A?—T 

For B2-M, f M , For M-B2, f = (1-19) 

In conjunction to the previous works by Liang (13), Rogers (16) and Tanaka (10), (21) presumed D,Q, H to be 
constant from the equation (1-12 and revised energy equation 

is given by 

G(T,a,0(. T ,a)eL 1 =G(» + G(7\0 (1-20) 

Later, a constitutive model was developed (22) on the basis of the principle of maximum transformation 
dissipation which is defined as 

D f (gY) = max CaeCaii D^C,-\), (1-21) 


Subsequently, the importance of elastic constants, a micro plane model is developed by (23) to configure a 
constitutive relation between the macroscopic and microscopic scale for the phenomolgical models, which yields to 


a-ij ^niUj dQ + ^^^(nAj + nj6 ri jdn 

£ ij =Tjn £ N n i n j dn+ hfn EJ T C n i S rj + n j S ri)dn 


( 1 - 22 ) 

(1-23) 


In compliance to the work by (23) a continuum model is proposed in (24) to model energy storage during 
martensitic phase transition, in terms of introducing internal variables such as martensitic volume fraction into evolution 
equations given by 

xp e = xp e (9, E e , z), ip s = xp s (9, Y ei j (1-24) 


From the above models, most of the work is based on uniaxial loading which do not have framework that can 
include orientation mechanism. Since it is difficult to predict the behavior of microscopic behavior based on macroscopic 
model, intact provision to study this phenomenon of microscopic level is by subjecting the specimen one dimensional, 
multiaxial proportional and non-proportional loading. An imperious work is established by (25), the specific free energy 
and dissipation equations for one dimensional, multiaxial proportional and non-proportional model are given by 

Uniaxial: pxp = pxp e (£,£ tr j + pxp T (Tj + pxp tr (£ tr ,T,Hj 

Multiaxial: pxp{§, T,£ tr ,H) —J§_— £ tr )- k: (e — £_ tr ) + pC |(T — T 0 ) — Tln (^jj + y[piz(T — T*j + H] (1-25) 


A plastic yield constitutive model is developed by (26) based on plastic yielding, which can estimate the 
simultaneously recoverable transformation strains and irrecoverable plastic strains. The Gibbs energy for the overall SMA 
material is given by 

G{p.T,E t ,S,g t ,EV,g v A ,gl l ,p) = (l - $jG A (<j,T,£r>,g v A ,p) + SG A {o.T,E*,g p w fi)+ G mix (a, £*,/) (1-26) 


An understanding from the literature can be observed that constraining of parameters is required to unify the 
thermo-elastic-constitutive equations by coupling at micro and macro-mechanics. After careful validations of the above 
models quantitatively by (27), (28), (29) and (30) etc., it is realized that different applications and characterization 
techniques yield different anomalies to the above models. Thus, new constitutive models have been adopted based on the 
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crystallography planes and their orientations. An outlook provided by (31) which decomposes the system into elastic part 
which represents stretching and rotation of overall structures and inelastic part, which is governed by generation, growth 
and annihilation of B2-M fine structure. In such case, the free energy is given by 

f(E e ,-9,0 = i£ e .C(f)[E e ] -09- <9 o M0f)C0f)[£ e ] + J09 - i9rX + - ctf - <9 0 ) - c<9 ln(J) (1-27) 

The paper also concluded the effect of strain rate influencing strain-hardening that it is due effect of thermal 
effects associated with phase transformations. A qualitative and quantitative correlation is made to understand the above 
concept. Extending the concept of (31), a constitutive equation for reorientation and detwinning phenomena is developed 
in (32) by mapping reorientation as movement of interfaces between habit plane variants (hpv’s) and detwinning as 
movement of interfaces between lattice correspondent variants (lcv’s). Neglecting the M-M interaction in the equation 
(1-27, the driving force on detwinning system and the driving force on hpv reorientation is given by 

t l = (C e T*).S l 0 , a ij = (C e r).Z‘ ; ' (1-28) 

An insight to consider the thermal elastic strain and inelastic strain due to accommodation of twins between 
martensite variants addition to (12) has been modelled by (33) and total macroscopic strain is given. 

£ = £ e + £ th + r + £ twin 


s e — S.a 

£ th = a{T - T ref ) 

E T - f£ T 

gtwin _ ^FAgtwin 

CONCLUSIONS 


0-29) 


From the above review, it can be presumed that shape memory effect (SME) is modelled microscopically, in a 
more precise way than the macroscopic hysteresis. The understanding of these alloys in the present paper restricts only to 
loading parameters to understand the hysteresis, but a more insight of exhibiting SME will be emphasized in the next paper 
by involving processing parameters. 
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